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Abstract 

We discuss the construction of wave packets resulting from the solutions of a class of Wheeler- De Witt 
equations in Robertson- Walker type cosmologies, for arbitrary curvature. We show that there always exists 
a "canonical initial slope" for a given initial wave function, which optimizes some desirable properties of 
the resulting wave packet, most importantly good classical-quantum correspondence. This can be properly 
denoted as a canonical wave packet. We introduce a general method for finding these canonical initial slopes 
which is generalization of our earlier work. 

PACS numbers: 04. 20. Ex, 04.60. Kz, 04.25.-g 

1 Introduction 

The question of the construction and interpretation of wave packets in quantum cosmology and its connection 
with classical cosmology has been attracting much attention in recent years. Moreover, numerous studies have 
been done to obtain a quantum theory for gravity and to understand its connection with classical physics. The 
problem of the relation between quantum cosmology and classical physics is an important one that exists even in 
simple models pQ. Most authors consider semi-classical approximations to the Wheeler-DeWitt (WD) equation 
and refer to regions in configuration space where these solutions are oscillatory or exponentially decaying as 
representing classically allowed or forbidden regions, respectively. These regions are mainly determined by the 
initial conditions for the wave function. Two popular proposals for the initial conditions are the no boundary 
proposal [2] and the spontaneous nucleation from nothing [3]. These proposals have been attractive to many 
authors because they lead to some classes of classical solutions represented by certain trajectories which posses 
important features such as predicting an inflationary phase. 

In quantum cosmology, in analogy with ordinary quantum mechanics, one is generally concerned with the 
construction of wave functions by the superposition of the 'energy eigenstates' which would peak around the 
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classical trajectories in configuration space, whenever such classical-quantum correspondence is mandated by 
the nature of the problem. However, contrary to ordinary quantum mechanics, a parameter describing time 
is absent in quantum cosmology. Therefore, the initial conditions would have to be expressed in terms of an 
intrinsic time parameter, which in the case of the WD equation could be taken as the local scale factor for the 
three geometry [J. Also, since the sign of the kinetic term for the scale factor is negative, a formulation of the 
Cauchy problem for the WD equation is possible. The existence of such a sign is one of the exclusive features 
of gravity with many other interesting implications. 

The construction of wave packets resulting from the solutions of the WD equation has been a feature 
common to most of the recent research work in quantum cosmology [5j [6] [7] . In particular, in reference [5] 
the construction of wave packets in a Friedmann universe is presented in detail and appropriate boundary 
conditions are motivated. Generally speaking one of the aims of these investigations has been to find wave 
packets whose probability distributions coincide with the classical paths obtained in classical cosmology. In 
these works, the authors usually consider model theories in which a self interacting scalar field is coupled to 
gravity in a Robertson- Walker type universe with zero k. The resulting WD equation is often in the form of an 
isotropic oscillator-ghost-oscillator and is separable in the configuration space variables. The general solution 
can thus be written as a sum over the product of simple harmonic oscillator wave functions with different 
frequencies. As usual, the coefficients in the sum and hence the exact form of the wave packets are determined 
by the initial conditions, which are usually specified by giving the wave function and its derivative at R = 0, 
where R is the scale factor. The coefficient are chosen such that the following desirable properties are obtained. 
Firstly, there should be a good classical-quantum correspondence, which means the wave packet should have 
compact support centered around the classical path, the crest of the wave packet should follow as closely as 
possible the classical path, and to each distinct classical path there should correspond a wave packet with the 
above properties. The first part of this condition implies that the initial wave function should consist of a few 
localized pieces. Secondly, one expects the square of the wave packet describing a physical system to posses a 
certain degree of smoothness. 

At the classical level, if the initial position is known, the corresponding momentum can be calculated from 
the constraint H = 0. However, in quantum cosmology the WD equation can be in general put in the form of 

2 



a second order hyperbolic PDE. Therefore, complete description of the initial condition would in general entail 
the specification of both the initial wave function and its first derivative. However, it is in general an unresolved 
problem as to how, by knowing the initial classical values of the positions and momenta, one may embark on 
obtaining the corresponding quantum distributions. Generally speaking, there is no generic way of doing this. 
What is sometimes done is to consider a distribution which peaks around the classical values of the variables of 
interest, mostly and for simplicity taken as a Gaussian distribution, and simply let the resulting wave packet 
evolve and hope for a good classical quantum correspondence (8j [9] . Indeed, even if the wave packet does 
follow the classical path it might not be smooth. Smooth wave packets have been constructed in the past by 
adapting coherent states for the initial wave function and the undetermined coefficient related to the initial 
first derivative distribution are chosen to have the same functional form as those of the initial wave function 
[10] . Here, we shall generalize the aforementioned prescription to be applicable to any reasonable potential 
in the same class of problems. Moreover, the prescription is general enough to be applicable to problems 
which are not exactly solvable. The situation is therefore becomes analogous to the classical case: Given an 
initial position distribution, there is a canonical initial momentum distribution which ensures good classical- 
quantum correspondence. This is in contrast to the usual practice where there is no uniform prescription for 
choosing initial momentum distribution in order to obtain a good classical quantum correspondence as possible 
is obtained. Our solution here is based on [lOj where a certain prescription was offered for obtaining initial 
canonical slope. 

The paper is organized as follows: In Sec. [5] we outline the main problem which is a case of Robertson- 
Walker type cosmology where the matter is taken to be a particular type of self-interacting scalar field. We 
derive the main equations both for the classical cosmology and the quantum one. In Sec.[3l we present a general 
prescription for finding a Canonical Initial Slope (CIS) for a given choice of suitable initial condition on the 
wave function, which produces the above desired properties of the wave packet. In this derivation we will come 
up with a decisive criteria for what initial wave functions are suitable for classical-quantum correspondence. 
This prescription is general enough to be applicable to both exactly solvable cases and the ones which are 
not. Therefore, we begin Sec. [3] with a description the Spectral Method (SM) [11] which is a robust numerical 
method. We then present the k = case which is exactly solvable and display the resulting wave packet 
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obtained by CIS. Then we consider the case k 7^ 0, for which neither the classical nor the quantum equations of 
motion are exactly solvable. In particular, in our model the classical equations of motion turn out to be a set 
of non-linear coupled ODEs with moving singularities. For the classical cases, we solve them by a numerical 
method that we introduced in an earlier work |12j . In the quantum cosmology case we solve the problem using 
SM. For this case we display the resulting wave packets obtained by the CIS prescription. In Sec. 2] we draw 
some final conclusions. 

2 Robertson- Walker cosmology with a coupled scalar field 

In Robertson- Walker cosmology, one often considers the coupling of a scalar field to gravity. The resulting 
field equations include a 'zero energy constraint'. The WD equation in quantum cosmology is the result of the 
quantization consistent with this constraint. Now we discuss our model and derive the field equations in the 
classical and quantum domains. 

Consider the Einstein's field equation coupled to a scalar field 

R vv - = KT^cf), (1) 

AV-f = 0, (2, 

where Tn U is the energy-momentum tensor coupling the scalar field <j> to gravity and U is a scalar potential for 
the scalar field <j>. Hereafter, we shall work in a system of units in which k = 1. We parameterize the metric as 

where R(t) is the usual scale factor and k — +1, 0,-1 corresponds to a closed, flat or open universe, respectively. 
The field equations resulting from {!]) and (|2|) with the metric given by (J3j) can be written as 
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Here a dot represents differentiation with respect to time. We require the potential U(<j)) to have natural 
characteristics for small 0, so that we may identify the coefficient of \<\? in its Taylor expansion as a positive 



mass squared m 2 , and U(0) as a cosmological constant A. An interesting choice of U(cf>) with three free 
parameters is © Q31 HD1 Q31 US] , 

b 



U(4>) = A 



— sinh (ac, 



2a 2 



smh(2ac! 



(7) 



In the above expression A may be identified with a cosmological constant, m 2 = d 2 U j dc\?\^ = Q is a mass squared 
parameter, b is a coupling constant, and a 2 = |. However, one of our main reasons for choosing this potential 
is that its exact solutions exist for both the classical and quantum cosmology for k = [lOj and we want to 
solve the problem for k ^ and compare the solutions to the corresponding ones for k — 0. To do this, we have 
general numerical methods that can solve both the classical [12] and quantum cosmological [11] cases. More 
importantly, having the exact solution for k = and approximate solution for k ^ for the same potential 
simultaneously, we can formulate a general procedure for finding the CIS and test it. Our claim would then 
be the resulting general formulation works for any potential in the same category. This gives us motivation to 
solve this problem for k ^= for the case of quantum cosmology. The Lagrangian giving the above equations 
of motion can be written as 



-3RR 2 + 3kR + R 6 [cp 2 /2 - U(4>)} 



(8) 



This Lagrangian can be cast into a simple form by using the transformations X = R 3 ' 2 cosh(a^) and Y = 
i? 3 / 2 sinh(a<^), which transform the term R 3 U(4>) into a quadratic form. Upon using a second transformation 
to eliminate the coupling term in the quadratic form, we arrive at new variables u and v, which are linear 
combinations of X and Y, 



cosh(6>) 
- sinh((9) 



- sinh(0) 
cosh(6») 



X 
Y 



where, 



1 ,_i/-26 
= - tanh 1 — tr- 



ill the terms of the new variables, the Lagrangian takes on the following simple form, 



L(u,v) = - 



1-1 1 2\ 1-2 2 2\ J i / 2 2\l/3 

(U — 0J x U ) — (V — ld 2 v ) — ~^k\ u — V ) ' 



(9) 



(10) 



(11) 
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where w 2 2 = — 3A/4 + m 2 /2 =p V m 4 — 4b 2 /2. The resulting classical field equations are, 

2 3/b IX - , 

u + u lU + — (u2 _ v2)2/3 = 0, (12) 

g + "fr + T ( u 2-V 3 = °' (13) 

,i 2 + w?u 2 - v 2 - lo 2 v 2 + h(u 2 - v 2 ) 1 '* = 0. (14) 

Equations (|12[) and (|13p are the dynamical equations and (|14[) is the zero energy constraint. The non-linearity 
and the moving singular behavior of these equations are now apparent. 

The corresponding quantum cosmology is described by the Wheeler-DeWitt (WD) equation written as, 

H4>(u,v) = j-J^ + + lo 2 u 2 - oj 2 v 2 + h(u 2 - v 2 ) 1 ^ i,{u, v) = 0, (15) 

which arises from the zero energy condition Eq. (|14p . 

3 Solutions for the classical and quantum cosmology cases 

In this section we shall present a general prescription for finding CIS. The cases that we shall study will include 
some cases which are exactly solvable (k — 0) and some which are not (k ^ 0). For the latter we have to resort 
in part to some numerical methods. Therefore, we start this section by a discussion of the numerical method 
that we shall use, so that later on we can concentrate on the physics of the problem with no interruption. The 
general PDE that we want to solve is 

where f'(u,v) is an arbitrary function. This equation has the general form of a hyperbolic PDE equation and 
is encountered in many branches of science and in particular physics. It is notable that such equations may 
represent a wave-like equation whose solution may rapidly oscillate. In such cases, the usual spatial integration 
routines such as Finite Difference Methods fail to produce a reasonable solution. It is therefore of prime 
importance to introduce a method of solution which is efficient, accurate, predictable and easily implemented, 
given any reasonable boundary conditions [11] . 



Equation (I15|) in general is not separable except when k = [10j . and we have to resort to a numerical 
method to solve this hyperbolic PDE. Among the various numerical techniques, we choose SM which has the 
following advantages: It is very simple, fast, accurate, robust and stable. Spectral Method (SM) [TBI [T71 [TB"] . 
consists of first choosing a complete orthonormal set of eigenstates of a, preferably relevant, hermitian operator 
to be used as a suitable basis for our solution. For this numerical method we obviously cannot choose the whole 
set of the complete basis, as these are usually infinite. Therefore we make the approximation of representing 
the solution by only a finite superposition of the basis functions. By substituting this approximate solution 
into the differential equation, a matrix equation is obtained. The expansion coefficients of these approximate 
solutions could be determined by cigcnfunctions of this matrix. In this method, the accuracy of the solution 
is increased by choosing a larger set of basis functions. Having resorted to a numerical method, it is worth 
setting up a more general problem defined by the following WD equation, 

d 2 d 2 



m(u,v) = ^-^ + -^+Lu 2 1 u 2 -Lu 2 2 v 2 + f(u,v)jT/j(u,v)=Q. (16) 

As mentioned before, any complete orthonormal set can be used. In this section we use the Fourier series basis 
by restricting the configuration space to a finite square region of sides 2L. This means that we can expand the 
solution as, 

2 

i>{u,v) = E ^2 A m,n,i,j 9 l (-J-) 9j » ( 17 ) 



where, 



91 M = ^ M ' and R m = ( ' (18) 
32 ( 2r F)=^Scos(^). I 2 > m = 

By referring to the WD equation (|16[) . we realize that in the Fourier basis it is appropriate to introduce /' as, 

f'(u,v) = f{u,v)+Lo\u 2 ~lo 2 2 v 2 . (19) 
We can make the following expansion, 

. , s , v-^ , /7717TU\ /niTV\ 

f( U , v)i>{u, «) = E E B m,n,i tj 5, ( — ) 9j ( — ) ■ (20) 
where B' m n i j are coefficients that can be determined once f'(u, v) is specified. By substituting (|17l20p in (fT6|) . 
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and using independence of gi (^ 2 p i )s and gj ( Il j^-)s we obtain, 



■^-m,n,i ,j + B m ni j 0, 



(21) 



where, 



B' 



E 



m ,n ,z ,,7 



L 
-L 



/mnu\ /nirv\ } , 
9i ) 9j {— ) J [u, v)g v 



m itu \ I n'mj 

gj> ( — I dudv 



A , , ; ■/ 



/ j ^m,n,i,j,m' ,n' ,i' ,j' -^m' ,n' ,i' ,y ■ 



(22) 



Therefore we can rewrite (1211) as 



/ ran \ 2 / wk \ 2 



A m ,n,i,j + ^ f C"m,n,i,j,m' ,ri ,i> ,j' "^r. 
m 1 7 n' ,i' j' 



0. 



(23) 



V L / \ L / 

In this case, we select 4iV 2 basis functions, that is m and n run from 1 to N. It is obvious that the presence 
of the operator f'(u, v) leads to nonzero coefficients C' m ni ■ m , n , j( ./ in ([23]) . which in principle could couple 
all of the matrix elements of A. Then we replace the square matrix A with a column vector A' with (27V) 2 
elements, so that any element of A corresponds to one element of A' . This transforms (f2"3")l to 



DA' = 0. 



(24) 



Matrix D is a square matrix now with (27V) 2 x (27V) 2 elements which can be easily obtained from (|2"3")) . Equation 
([24| can be looked upon as an eigenvalue equation, i.e. DA' a = aA' a and the matrix D has (2N) 2 eigenvectors. 
However, for constructing the acceptable wavefunctions, i.e. the ones satisfying the WD equation (|16p . we only 
require eigenvectors which span the null space of the matrix D. That is, due to (|23|) we will have exactly 27V null 
eigenvectors which will be linear combination of our original eigenfunctions introduced in (|17[) . After finding 
the 27V eigenvectors of D with zero eig envalue, i.e. A' k (k = 1, 2, 3, 27V), we can find the corresponding 



elements of matrix A, A^ ni -. Therefore, the wavefunction can be expanded as 



4>(u,v)=J2^ k *P k (u,v)=J2 Xk J2 A kn,i,j 



U7TV 



L J \ L 



9j 



(25) 



Here A* s in (|25[) are arbitrary complex constants to be determined by the initial conditions. 

We are free to adjust two parameters: 27V, the number of basis elements used and the length of the spatial 
region, 2L. This length should be preferably larger than spatial spreading of all the sought after wave functions. 



However, if 2L is chosen to be too large we loose overall accuracy. Therefore it is important to note that for 
each TV, L can be properly adjusted [16] , 

The most important physical aspects of this problem that we want to address is on the question of the 
existence of a CIS. In our previous investigation of the same problem for k — [10] we discovered that a 
"canonical" choice for the slope exists whose use produces wave packets with all the desired properties stated 
in the introduction. Here we have come up with a more general and systematic method for obtaining this 
slope. The results are based on the same general prescription as before and obviously give identical results 
for comparable cases. This slope turns out to depend on the nature of the hamiltonian and the initial wave 
function. Since the WD equation is a hyperbolic PDE, for any reasonable choice of initial wave packet, 
there will remain certain undetermined coefficients which will be pending the specification of the initial first 
derivative. Our prescription is simply to choose the undetermined coefficient to have the same functional form 
as the determined ones. We shall clarify this prescription further in some concrete examples. If the problem 
is exactly solvable this slope can be easily obtained by the method described in Ref. [10]. If it is not, we can 
use the following prescription: We use the minimum approximation necessary to make the general PDE (i.e. 
(fT6|) ) separable near v = and solve the resulting equations. It is obvious that the presence of the odd terms 
of v dose not have any effect on the form of the initial wave function but they are responsible for the slope of 
the wave function at v = 0, and vice versa for the even terms. We can approximate (| 1 5|) near the v = 0, so up 
to the first order in v we have, 

()1 + ^- +UJ lu 2 + -Mu 2 ) 1 lAi,{u,v) = Q. (26) 



du 2 dv 2 4 
This PDE is separable in u and v variables, so we can write 



tp(u,v) = (p(u)x{v). (27) 



By using this definition in ([2"6"jl . two ODEs can be derived, 

d 2 Xn(v) 



, , . E nXn (v) = 0, (28) 
civ 



d 2 ip n (u) ( o 2 9 2/l\ I \ „ , , 

- — ^2 VU^lu + -ku 1 J tp n (u) = E n ip n (u), 



(29) 



where E n s are separation constants. These equations are Schrodinger-like equations with E n s as their 'energy' 



levels. Equation (128)) is exactly solvable with plane wave solutions or, 

X„(v) = a n cos {^/E^ v^j + i/3 n sin (JeZ v^j , (30) 

where a n and [3 n are arbitrary complex numbers. Equation (|29p does not seem to be exactly solvable and we 
resort to a numerical technique. As mentioned before, SM can be used to find the bound state energy levels 
(E n ) and the corresponding wave functions (ip n (u)) with high accuracy. The general solution to the (|26p can 
be written as, 

tp(u,v)= ^ {A n cos(\jE^v)+iB n sin{y/E^v))tp n (u)+ ^ (C n cos(y^v)+iD n sm(^/E~v))ip n (u). (31) 

ra=even n=odd 

The separation of this solution to even and odd terms, though in principle unnecessary, is crucial for our 
prescription for the CIS. As stated before, this solution is valid only for small v. The general initial conditions 
can now be written as, 

n + J2CnVn{u) (32) 

even odd 

tp'(u,Q) = i ^ B n y/E^ cp n (u) + i^D n \[~E~ n Lp n (u) , (33) 

even odd 

where prime denotes the derivative with respect to v. Obviously a complete description of the problem would 
include the specification of both these quantities. However, given only the initial condition on the wave 
function, we claim there is only a CIS which produces a canonical wave packet with all the aforementioned 
desired properties. We can qualitatively describe our prescription for this case as setting the functional form 
of the odd undetermined coefficients to be the same as the even determined ones and vice versa. This means 
that the coefficients that determine CIS i.e. B n for n even and D n for n odd, are chosen by the following, 

B n = C n for n even D n = A n for n odd (34) 

To be more specific, although C n (A n ) is defined only for n odd (even), we can extend its definition to n even 
(odd) by choosing the same functional form. 

Now we want to discuss the general settings for quantum-classical correspondence. As mentioned before 
in the classical cosmology case, the initial condition is given by specifying a few isolated initial points in the 
configuration space and uniquely determined velocities associated with each one through the constraint H = 0. 
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In order to established a good quantum-classical correspondence, the first thing we should do is to specify 
an initial wave function with isolated peaks centered around the classical ones. Then we can use the CIS 
prescription to obtain the initial first derivative distribution. One of the main outcomes of our prescription is 
that the expectation value of initial momentum computed using our prescription in the quantum cosmology 
case, denoted by [(p v )] v= qi is approximately the same as the corresponding classical quantity [i)(f)] (=0 (in our 
notation) for each isolated classical point. That is [(Pt>)] v _ — ¥ [v{t)] t=0 where, 



Obviously the domain of the above integral should extend only over the appropriate part of the initial wave 
function, which corresponds to the classical point in question. This relationship approaches an equality for 
initial states for which classical description is appropriate. We shall shortly make this statement more explicit 
for the particular problem under investigation in this paper. Our prescription for the CIS is general enough to 
be applicable to any initial wave function. 

In the case k — the problem is exactly solvable [10] and has a closed form solution in the Simple Harmonic 
Oscillator (SHO) basis. In references 0[TO] the authors considered the following initial condition for the wave 
function, 



which consists of two symmetric coherent states of a one dimensional SHO. Note that, for the given choice of 
the initial wave function we can not expect a good classical-quantum correspondence for % < 3, since there 
would be a significant overlap between the two pieces, which could cause quantum interference. With this 
choice for the initial wave function one might expect the classical-quantum correspondence to be manifest for 
X > 3. However, as we shall see, this crucially depends on the choice of the initial slope. Obviously, this choice 
for tp{u,0) is not unique and is chosen so that the resulting wave packet would have compact support. In the 
case k = 0, (|15|) is exactly separable without any approximation and yields exactly the same Schrodinger-like 
equation, i.e. SHO, as the one obtained from (|29p for both u and v variables with frequencies lo\ and u>2 
respectively. The solutions to these equation are well known, i.e. for u is 




(35) 




(36) 






(37) 
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where H n s are the Hermite polynomials and superscript zero indicates k = 0. Then the complete eigenstates 
will be products of eigenstates for u and v, with constraint that they should have the same energy. In this 
case the initial wave function is chosen to be even and taking into account our prescription for the CIS (|34p , 
(|32I33[) can be rewritten as, 

4>(u,0) = 5>„^°( W ) (38) 

even 

^'(u,0) = i£W(2n+l)u;i (39) 

odd 

where the choice, 



= P -k\x\*J£ 



(40) 



V2"d ' 

produces our aforementioned wave function (|36[) and x is an arbitrary complex number. The CIS introduced 
in Ref. [10], based on the same prescription as introduced here but for the exactly solvable case k = 0, is, 

i/> («, 0) = 2^ ( -i)(~/2) n! . ( 41 ) 

"odd (n/2)! 



where prime denotes differentiation with respect to v. Numerical comparison between (|39|) and (|4Tj) shows that 
their relative difference is O (10 -4 ). Their difference is due to the fact that (|4"Tj) is exact and is approximate. 
We can now easily construct the coherent wave packets using (j3"B"|) and ([31?]) or (|4"Tj) . For illustrative purposes 
we choose the simplest case possible which is when toi = u>2- The result is shown in the left part of Fig. [TJ As 
can be seen from this figure the parameters of the problem are chosen such that the initial state consists of two 
well separated peaks and this class of problems are the ones which arc also amenable to a classical description. 
We will make this statement more quantitative at the end of this section. We should mention that there are a 
variety of different cases illustrated in Ref. [TU] for the same exactly solvable case of k = 0, including u>\ ^ ll>2 
obviously using (|4"Tj) . 

The classical paths corresponding to these solutions can be easily obtained from (|12I13[) . The corresponding 
initial conditions for the classical case is the following: 

«(0) = -X. «(0)=0, «(0) = 0, 6(0) = 6o, (42) 

and the parameters x an d are adjusted so that the zero energy condition (fT4"|) is satisfied. With this choice 
of initial conditions the exact classical paths would be, in the case k = 0, Lissajous figures in general, which 
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are circles with radii x when lo\ = L02 = u>. Using (|35() we can compute the average initial momentum in the 
quantum case for the left peak, where the corresponding classical values are u(0) = — x- w (0) = 0, u(0) = 
and v = +WX- It is important to note that the domain of integration in equation (135[) in this case should be 
from — 00 to zero, as explained before. We obtain the following two expressions first using the exact value for 
CIS pTj) . and second the approximate CIS 



r, xiexact . \p C n C m yn H ^;:yujn mV J ) ^ a n \v ) m +n-2k-V 



(m/2)!Vw^(0) 



(-1) ((m+l) /2) m J 



Min{m,n} , 



2 "i?(0) m +„_2fe- 



fc!(n — fc)!(m — fc)! 



„ „ r , Mi„{m,n} ^krrf^ 

\(n \l a PP rox - - 4 V V iJ(2m+l)uj n'm' V f^Jy+n-gfe-i ...x 

Having precisely set the initial conditions for both the classical and quantum cosmology cases, we can now 
superimpose the results for the case fc = 0, as illustrated in the right part of the Fig. [1] As can be seen 
in the figure, and also for all the cases presented in [10] . the classical-quantum correspondence is manifest. 
In Fig. [2] we show v(Q) versus x f° r the three cases: classical, exact CIS, and approximate CIS. As can 
be seen from this figure the classical-quantum correspondence is meaningful only for % > 3, that is where 
there is no significant overlap between the two pieces of ip(u, 0). This correspondence becomes better as x 
increases. That is [{Pv)] v= q = [^(^)]( = o = ^0 for large x- We can generically expect loss of classical-quantum 
correspondence when there is significant quantum interference between different segments of the wave packet 
in the configuration space, when each segment is supposed to correspond to a distinct classical configuration. 
We would like to emphasize that the choice of CIS is crucial for establishing a good correlation. For example 
in Ref. [5 the initial slope in the quantum case is chosen to be zero, implying zero initial average momentum 
(see (|35jl ). however the corresponding classical quantity v(0) — ±cux an d the classical path would again 
be a circle. To be more precise, the WD equation emanates from the zero energy condition and their choice 
of initial conditions violates this constraint. Therefore it would be impossible to establish a classical-quantum 
correspondence in this case. It is interesting to note that their resulting wave packet (\ip(u, v)\ 2 ) is wildly 
oscillating. 

For the case fc ^ 0, the problem is not exactly solvable in quantum cosmology and we will use the SM to 
get an approximate solution. We have to mention that the corresponding equations for classical cosmology 
(fT2)) - (ri"4")) . are a set of nonlinear, coupled ODEs with moving singularities which are not exactly solvable either. 
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However, a general method for solving them has been presented in [12j . and this is the method we shall use. Also 
a detailed explanation of the physical setting of the problem in the classical domain has been presented in [19] . 
For ease of comparison, we choose the same illustrative problem with u>i = lo-2 = lu and the same coefficients for 
the initial wave functions as the case k = (i.e. (|38l40p ). However, note that since the eigenstates and 
4>n( u ) (both obtainable from (HU)) are slightly different, the resulting initial wave functions for k = ±1 would 
be slightly different from (|36| . From (|29I33[) it is apparent that the CISs of the wave functions are different in 
all cases (k = 0, 1, —1), due to the differences both in 4>n(u) s and E n s. We have computed both the classical 
and quantum solutions for the k = +1, —1 cases and the results are shown in Figs. [3] and 01 respectively. Note 
that the wave packets are very smooth and we have good classical-quantum correspondence for both cases. 
Moreover, we have computed the graphs of v(0) versus \ f° r & — an d have found the classical-quantum 
correspondence to be as good as that exhibited in Fig. [2] for k — 0. In order to highlight the differences 
between the classical paths (and obviously the wave packets) we exhibit r = \/u 2 + v 2 versus 9 = arctan(ti/u) 
in Fig. [5] for the three cases. 

4 Conclusions 

We have described a Robertson- Walker type cosmology leading to classical dynamical equations given by ([T2"l 
03}, and the corresponding WD equation represented by ([T5|) . All these equations are exactly solvable for the 
particular potential chosen when k = [10], and do not seem to be so when k ^ 0. In the latter cases we have 
solved these equations numerically: the classical cosmology equations by the numerical method introduced 
in [12j . and the quantum cosmology case by an implementation of the SM [11) . Most importantly, we have 
introduced a general method for finding the canonical initial slopes for any initial wave function and any value of 
k, whose use produces wave packets which can be named canonical wave packets. The latter have the following 
desirable properties: The correlation part of classical-quantum correspondence is acceptable, as outlined in 
detail in the Introduction. Moreover, the resulting wave packets are very smooth and this is usually considered 
to be a desirable property. Our investigation has also revealed that the classical-quantum correspondence is 
possible only when there is no overlap between the lumps of a suitably chosen initial wave function. Using the 
CIS, the average initial momenta become very close to their corresponding classical values in the limit where 
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the overlap between the separate parts of the initial wave function goes to zero. It is important to note that 
our generalized method is also applicable to cases which are not exactly solvable. 
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Figure 1: k = case: Left, the square of the wave packet \ip(u, v)\ 2 for x = 3.5 and N = 15. Right, the contour 
plot of the same figure with the classical path superimposed as the thick solid line. 
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Figure 2: v(0) versus % for the three cases: classical (dashed-dot), exact CIS (dashed), and approximate CIS 
(solid). On this scale the latter two approximately coincide for \ > 4. Their maximum difference is of order 
1CT 2 . These two CISs asymptotically approach the classical value for \ > 8-) 



17 



2 4 6 



-6 -4 -2 
u 







Figure 3: k = 1 case: Left, the square of the wave packet \ip(u, v)\ 2 for x = 3.5 and N = 15. Right, the contour 
plot of the same figure with the classical path superimposed as the thick solid line. 
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